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Abstract. IRS 13E is an infrared, mm and X-ray source in the Galactic Centre. We present the first Chandra 
X-ray spectrum for IRS 13E and show that it is consistent with a luminous and highly absorbed X-ray binary 
system. Since the X-ray luminosity is too large for a solitary star, our interpretation is that of an early-type 
long-period binary with strong colliding winds emission. This naturally explains the observed X-ray spectrum and 
count rate as well as its lack of significant short term variability. Due to the short lifetime of any nebula 0.2 pc 
from the putative central super-massive black-hole, we argue that the primary of IRS 13E has exited the LBV 
phase in the last few thousand years. 

Key words. X-rays: stars - Stars: Wolf-Rayet - Stars: individual: IRS 13E 


1. Introduction 


It is probable that Sgr A*, the compact, nonthermal ra¬ 
dio source at the Galactic Centre (GC) is a 2-3x10^M(r ) 
black-hole (for a recent review see [Melia fc Falcke 2001 ). 
Pervading the central parsec of the Milky Way is a clus¬ 
ter of a few dozen Hel and early-type stars (^ellgren et al. 


|199C ; Genzel et al. 1996| ). One of these stars, IRS 13E, has 


been identified as a Wolf-Rayet (WR) star with spectral 


class WNIO (Najarro et al. 1997) and lies within IRS 13, 


a compact HII region. The IRS 13 complex, dominated by 
IRS 13E, has been identified as a Hel, Pa-a, [Felll], and 


Hell line source ( 

Lutz et al. 1993; 

Libonate et al. 1995; 

Krabbe et al. 199£|; Stolovy et al. 1999 

)■ 

Motivated by the conjecture ( 

Coker & Pittard 2000) 


that IRS 13E is an X-ray binary, we present here an anal¬ 
ysis of Chandra Advanced CCD Imaging Spectrometer 
(ACIS) GC observations (obs ID=242) of this source. 


Details of the observations can be found in Baganoff et al. 
|(2000| , ^OOla The X-ray source ~ 4" west-southwest of 
Sgr A* is the only source in the central parsec that has 
been associated with a previously known stellar object, 
IRS 13E. Since IRS 13E appears as an X-ray source while 
the other early-type stars in the central parsec do not, 
IRS 13E must harbor a distinctive object. 

Based on its lack of significant variability at all wave¬ 
lengths and its strong X-ray luminosity with characteris- 
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tic kT ~ 1.0 keV, we argue that IRS 13E is most likely 
an early-type wide binary system with the primary only 
recently having exited the luminous blue variable (LBV) 
phase of evolution. 


2. Source Identification 


As discussed in poker &: Pittard (2000] ), we identify the 
X-ray source 3" west^J and 1.5" south of Sgr A* with the 
infrared source IRS 13E. However, further identification 
is far from clear. In K-band observations the IRS 13E 


complex is resolved into 3 components (Ott et al. 1999): 


IRS 13E1 {niK = 10.26), IRS 13 E2 (my = 10.27), an d 
IRS 13E3 {rriK = 10.31) (Fig. 8 of Paumard et al. 2001 ). 
However, these objects do not match up with mm observa¬ 
tions, which show a different set of 3 component s, some of 
which ar e extended (Fig. 2 of Zhao fc Goss 199^ ) . We pre¬ 
viously ( Coker fc Pittard 2000 ) identified the mm sources 
IRS 13E and IRS 13W with the IR sources IRS 13E1 and 
IRS 13E2, respectively, but careful astrometry shows that 
none of the K-band sources are coincident with the mm 
sources (Maillard 2001; private communication). This im¬ 
plies the mm sources may not be stellar; the relative lo¬ 
cations of the mm and IR sources open the possibility 
that the radio sources are due to the colliding winds of 
a binary - or even triple - system. However, it appears 
that in agreement with Paumard et al. (2001|), IRS 13E3 


At a distance of 8.0 kpc (Reid 1993), V ~ 0.04 pc. 
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is slightly extended so its identification as a single stellar 
source is unclear. _ 

Paumard et al. (2001) identify IRS 13E3 as the Hel 
source while Majarro et al. (1997| ) identify IRS 13E1 as 
the Hel source. However, NICMOS data, taken in March 
1998, show (Stolovy 2001; private communication) that 
IRS 13E1 has very little (if any) excess emission at 1.87/im 
while IRS 13E3 has some excess and IRS 13E2 has a strong 
excess. Since a 1.87/im excess probably represents a blend 
of Pa-a, Hel and Hell emission lines, IRS 13E2 is very 
clearly the dominant Pa-a and Hel source of the IRS 13E 
complex. 

Paumard et al. (2001) also report that their June 
1998 observations show IRS 13E3 is slightly elongated and 
much dimmer (rriK = 11.73) than the other two compo¬ 
nents. But the series of observations by Ott et al. (19991) , 
undertaken from 1992 to 1998, show that all three com¬ 
ponents are of nearly equal magnitude and that IRS 13E3 
(their source 21) in particular is non-variable at the level 
0.2 mag. Further, K-band images taken in 1993 and 1994 
( pfamblyn et al. 1996 ) which do not resolve the IRS 13E 
complex have ttik = 9.1 mag, fully consistent with the 


3 equal sources with rriK = 10.3 mag seen by Ott et al. 
|(1999| ). The cause of this ~ 1.5 mag discrepancy is unclear. 

Although not critical to our results, in this paper we 
will identify the WNIO star with IRS 13E2 and the X- 
ray source with the IRS 13E complex as a whole. The 
components of IRS 13E are separated from each other 
by ~ 0.1-0.2" in projection, so that, given the ~ 0.5" 
resolution of Chandra, we lack sufficiently good positional 
references for precise identification of the X-ray source. 


3. The X-ray Data and Spectral Fitting 


Chandra observed the GC in September 1999 and again 
in October 2000 ( Baganoff et al. 2000| , 2001a). Here we 
present detailed analysis of only the first epoch of ob¬ 
servations. We used standard data processing tools avail¬ 
able from the Chandra X-ray Center (CXC) as part of 
the CIAO v2.1 package^ to extract spectra of the X-ray 
source and of adjacent background, and to determine the 
appropriate response matrix and ancillary response (effec¬ 
tive area) function for the source extraction region. The 
source spectrum was extracted by binning events lying 
within a 2.5" radius circle centered on IRS 13E; the back¬ 
ground spectrum was extracted from within three 2.5" 
radius circles adjacent to and immediately north, south, 
and west of this source region. We avoided the region im¬ 
mediately to the east of IRS 13E, as it is dominated by 
intense emission from Sgr A*. Within the source region 
we find 137 counts in 46 ksec net integration time. The 
background regions contain ~ 45 counts in an equivalent 
area, such that the background-subtracted count rate is 
0.002 cts s“^. Within the limited Poisson statistics, the 
X-ray flux from IRS 13E appears to remain constant be- 


For details see jittp:/ /cxc.harvard.edu/ciao/index.htm] 



Fig. 1. A plot of counts per second per keV versus energy 
for the data (crosses) and best-fitting model (solid line). 


tween the two epochs. In addition, within the two ~ 50 
ksec observations, IRS 13E is not seen to vary. 

We fit the data with the MEKAL (Mewe et al. 1985) 
model in XSPEC0. This is a thermal plasma model and 
assumes optically thin X-ray line and continuum emission. 
We account for an external absorption column by simul¬ 
taneous application of the wahs model. 

In fitting the spectral X-ray emission, we vary three 
parameters: the X-ray luminosity in the 0.2-10 keV band 
(Lx), a characteristic temperature (fcT) of the X-ray emit¬ 
ting gas, and the line-of-sight column density (A^h)- The 
metallicity, Z, of the stars in the central parsec is not 
well known; ISM gas phase abundances are roughly twice 
solar ( Simpson et al. 1995) while that of GC red super ¬ 
giants is only solar ( |Carr et al. 200C ; Ramirez et al. 2000 ). 
The goodness of fit does not change significantly when 
the global abundance is varied so we have assumed solar 
metallicity for simplicity. 

Spectral models can b e htted to d ata with low counts 
using the Cash statistic ( [Cash 197£ ). However, a back¬ 
ground cannot be subtracted using this method. We 
therefore employ the fit statistic to our background- 
subtracted data, rebinned to a minimum of 10 cts per 
bin. 


4. Results &. Discussion 

Fig. ID shows the ACIS data along with the best-fitting 
model, which yielded xt — 1-5, formally a somewhat 
poor fit. The best-fitting parameters are listed in Table 
Inspection of Fig. shows that our model underestimates 
the evident line emission at ~ 3 keV and slightly overes¬ 
timates the emission at ~ 3.75 keV. Some of this may be 
due to poor statistics but the strong emission near 3 keV 
may also reflect a complex metal abundance. 

Our ht suggests A^h — 15 x 10^^ cm“^. This is 

somewhat larger than the 5-10 x 10^^ cm“^ found for 
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Table 1. Computed Model for IRS 13E and Arches A2 



IRS 13E 

Arches A2“ 

z 

1.0 

2.6*’ 

L. (Lq) 

7 

20 

kT (keV) 

1.0 ±0.4 

1.0 ±0.5 

Ah (10^^ cm"^) 

15 ±5 

12 ±2 


a) from Yusef-Zadeh et al. (2001b) 

b) averaged Z over Si, S, Fe, Ar, and Ca 


Sgr A* ( Baganoff et al. 2000 , 2001a ), implying that there 
is additional substantial absorption close to the IRS 13E 
X-ray source. If so, and if the size of the X-ray source is 
comparable to the ~ 0.1” diameter region seen in mm 
observations, then the characteristic density of the gas 
surrounding IRS 13E is consistent with a colliding wind 
binary (CWB) system (Stevens et al. 1992) but some¬ 
what more dense than LBV ejecta ( Reward et al. 2001 ). 
However, some of our estimated column to the GC may 


be due to dust (Baganoff et al. 2001a) since, unlike most 
of the massive stars in the central parsec, the IRS 13 com¬ 
plex is enshrouded by warm dust ( pLieke 1999 ). Similar 
large amounts of dust are seen around “cocooned” stars 
located in the Quintuplet Cluste r that is 50 pc from Sgr A* 
in projection ( Figer et al. 1999| ); although these s tars are 
proba bly W Cd type stars, this is not yet certain ( |Monetj 
et al. 2j^. 

The fitted characteristic temperature, kT ~ 1.0 keV, 
is consistent with the shocked winds of a CWB system 
such as a WR-bO. In contrast, an accretion source with 
a compact companion such as in a massive X-ray binary 
(MXRB) system, would typically have a harder spectrum 
(e.g., pchlegel et al. lO"^ ). Also, a solitary massive O- 
star typically has a characteristic temperature of only 
~ 0.5 keV ( Chlebowski et al. 1989| ). In the case of a CWB, 
the characteristic temperature represents a global average 
of the hot, shocked colliding winds, while X-ray emission 
from a MXRB probes the accretion disk. On the other 
hand, the X-ray emission from LB Vs is not well known; 
only rj Car, a suspected binary whose primary is an LBV, 
has a well-determined spectrum. The emission from t] Car 
is a combination of hard {kT ~ 5 keV) compact emission 
from the sta r and softer {kT ^ 0.5 keV) extended emission 
from ejecta ( Reward et al. 2001 ). However, as a binary, the 
X-ray emission from r] Car may not be typical of solitary 
LBV stars. 

The absorption-corrected X-ray luminosity of IRS I3E 
in the 0.2-10 keV band is found to be Lx — 7 Lq. Although 
larger luminosities are found for the brightest early-type 
binary systems (e.g. WR 140 and, assuming binarity, rj 
Car; van der Hucht et al. 1994; Corcoran et al. 2001), 
Lx — 7 Lq is still brighter than the typical CWB system. 
However, this is not wholly unexpected given t he large 
mass-l oss r ate and wind velocity of IRS I3E2 ( Najarro 
et al. 1 997 ). In contrast, such an X-ray luminosity is a bit 
low for a MXRB with an accreting compact source unless 
the system has an unusually long period. Additionally, 


the X-ray luminosity of a MXRB generally varies on short 
time-scales. Therefore, due to its low luminosity, lack of 
variability, and relatively soft spectrum, it is somewhat 
unlikely that IRS 13E contains a short period MXRB. 

The lack of variability also makes it unlikely that the 
X-ray emission from IRS 13E is due to either a flaring 
proto-star or young stellar object (YSO). The estimated 
intrinsic X-ray luminosity is also considerably higher than 
any known YSOs (e.g., Garmire et al. 2000). However, 
plenet et al. (200l] ) show that the K-L colour of the IRS 13 
complex as a whole is consistent with the presence of 
YSOs. 

Single massive stars seem to obey the rough rela tion 
logio (^x^^/^boi) = — 7 ± 1 ( Pallavicini et al. 1981 ), al- 
though more recent work ([Chlebowski fc Garmany 1991 


Moffat et al. 2001 ) suggests the scatter may be more sub¬ 
stantial and the ratio may be half a dex larger. L^®^ is 
the X-ray luminosity corrected for extinction due to the 
ISM but not for intrinsic extinction ( jWaldron et al. 1998 ). 
If we assume that = 10^^ cm“^, then our model for 

IRS 13E results in logio (Ax^^/^boi) = —7, no greater 
than solitary X-ray sources. CWB systems tend to have 
an enhanced X-ray to bolometric luminosity ratio by a 
factor of a few compared to solitary stars, but there is 
considerable variation. Possible contamination and large 
error bars in determining logio (Lx^’^/Lboi) for IRS 13E2 
make it difficult to draw any conclusions concerning bina¬ 
rity using this ratio. For example, it may be that very little 
of the extinction is intrinsic to the IRS 13E system (see 
below); this would result in a substantially larger L^®^. 

We must caution that given the low signal-to-noise of 
the data, the best-fitting parameters are not very well con¬ 
strained. For example, forcing kT = 5 keV and varying 
only A^h and Lx results in a fit with xt = 1-9- However, 
based on observations of other objects in the GC it is 
probable that ^ 10^^ cm“^. Given this additional 

constraint, we can say that Lx ILq and kT ^ 1.8 keV. 

Also shown in Table ^ are the results for the bright 
soft component of sou rce A2 of the Arches clus ter located 
~ 25 pc from Sgr A* ( Yusef-Zadeh et al. 2001a| ,[b|). The fit 
to this source implies a higher than solar metallicity for 
a range of metals. Although we assume solar metallicity 
for IRS 13E, the peaks near 3 and 4.25 keV in Fig. ^ are 
possible indicators of high S and Ar content. Enhanced 
abundances for these elements lead to a slightly better 
fit (Xy = 1.3) but due to the low number of counts in 
our spectrum, we do not attribute much significance to 
this. Our column density, characteristic temperature, and 
X-ray luminosity are close to that for A2, which is also 
coincident with a known IR stellar source, suggesting the 
objects are of similar nature. 

In short, even given the large uncertainties, the X- 
ray luminosity and characteristic temperature of IRS 13E 
do not favour a single object (0-star, WR star, LBV, or 
YSO). Also, the lack of variability is inconsistent with a 
short-period MXRB. On the other hand, all of the charac¬ 
teristics of IRS 13E are fully consistent with a long-period 
CWB system. 
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5. Evidence for a Recent Post-LBV System 

Although many significant details of massive star evolu¬ 
tion remain to be addressed, the present picture of mas¬ 
sive star evolution is that stars with zero-age main se¬ 
quence (ZAMS) mass ^ 25 Mq pass through the WR 
stage. The evolutionary sequence for the more massive 
stars (Mzams ~ 40 Mq) is thought to be (Walborn 1989; 
hanger et al. 199^ ): 


O ^ H-poor WN ^ LBV ^ H-free WN ^ WC ^ SN . 

Most stellar evolution models suggest it takes more than 
5 X 10® yrs for a massive star to evolve through the 
entire WR sequence, with the LBV phase taking more 
than 10'^ yrs ( ^tothers fc Chin 1996| ). While in the LBV 
phase, massive stars are thought to repeatedly move back 
and forth across the HR diagram, being bluer in quies¬ 
cence and redder during major eruptions. More than 70% 
of the LBV phase is spent in quiescence ([van Genderen 


2001). However, some evolutionary models suggest that 
extremely massive stars (Mzams ~ 60Mq), particularly 


those with high metallicity or large rotation rates (Meynel 
fc Maeder 2000 ), may neve r become true LBVs ( prowthCT 
et al. 1995 ) or even WRs ( Mowlavi et al. 1998 ). 


Since IRS 13E2, the presumed Hel line source (see §2) , 
is spectrally a late-type WN star ( van der Hucht 2001]) , 
it can be either a pre- or post-LBV object. As chemi¬ 
cally t he source ha,s a, large He to H ratio (Najarro et a,1. if is not significantly variable in (Zhao fc Gos; 


also have He/H of less than a few (Crowther et al. 1995), 
it seems unlikely, even in the peculiar environment of the 
GC, that 7 LBVs would exist at the same time in the same 
cluster. Comparison with WR 122, a star still shrouded in 
ejecta and thought to be just post-LBV, might be useful. 

Models of the bolometric luminosity and effective stel¬ 
lar temperature of IRS 13E2 (Mboi = —II.2 mag and 
Teff ~ 29 kK) place it on the left edge of the Humphreys- 
Davidson (HD) line, the location on the HR diagram 
where LBVs are thought to exist. This temperature, as 
well as the temperature of the other early-type stars in the 
central parsec, is somewhat colder than the lower limit of 
traditional WR stars (30kK; Wolf 198"^ ), but, due to the 
high extinction, determining effective temperature in the 
GC is notoriously difficult ( Dessart et al. 1999 ). At present 
only a few central stars have estimated Mboi and Teg val¬ 
ues; within the relatively large errors, all fall below or are 
near the HD line. The effect of metallicity on the HD line 
is unclear but models suggest that high metallicity moves 
LBVs to lower Lboi and Teg ( ^tothers fc Chin 1996 ). In 
the future, observations of CO absorption at 2.3/im and 
H 2 O absorption at 1.9/rm may better determine Teg and 
Mboi for IRS 13E2 and the other early-type stars in the 
GC (Blum et al. 1999). 

An LBV is often variable across many time-scales and 
frequencies. Although IRS 13E2 may be marginally var i- 
able at the level of 0.1 mag in K-band (Ott et al. 1999)^, 


T^), IlHS 1.3E2 is probably a post-LBV object. There' 1998D nor apparently in the X-ray ([Baganoff et al. 2001a| ). 


are thought to be a minimum of 6 WC^ stars in the cen¬ 
tral cluster ( van der Hucht 20'0l| ) so at least some cluster 
members, presumably those with the largest ZAMS mass 
if the cluster is coeval, have already gone through the LBV 
phase. This is particularly noteworthy since the ZAMS 
mass of IRS 13E2 is thought to be as high as ~ 120 Mq 
( ^challer et al. 1992 ), implying that the ZAMS mass of 
the 6 WC stars is even higher. Some of the Hel stars in 
the central cluster are comparatively high in H; they may 
be pre-cursor LBVs. If the statistics of LBVs versus WRs 
in the LMC is comparable to that in the GC, one would 
expec t 1 LBV out of a population of a few dozen WR 


stars ( [lmn£hre2^_&_^avldsonM994Parke£2997 )^ 

The H-poor WN stars, being pre-LBV, are more mas¬ 
sive, more luminous and larger than the H-free WN stars. 
As the terminal velocity of radiatively driven winds from 
early-type stars is in general correlated with their surface 
gravity, the H-poor WN stars should have slower winds 
and narrower lines than the H-free WN stars (see, e.g.. 
Smith et al. 1996; Conti 19^). Thus, the 7 narrow-line 


This suggests it has ended its LBV phase, although it may 
also be merely quiescent. Since it takes ~ 3 x 10® years 
for a 100 Mq star to proceed from the Z AMS through 
to the final WR phase ( Langer et al. 1994 ), this gives an 
approximate age of the central cluster. 

The high terminal wind velocity (cs 1000 km s“^) and 
lack of short-term variability suggests that IRS 13E2 has 
ended its LBV phase. Tidal disruption by Sgr A* means 
that any nebula surrounding IRS 13E2 is likely to be short 
lived after the star becomes a broad-line H-free WN star. 
In addition, since LBV terminal wind velocities can be 
lower than ~ 200 km s“^ (Najarro et al. 1997; Parker 


1997), any GC LBV-spawned nebula will be distorted by 


stars that Paumard et al. (2001) classify as LBVs could 
be H-poor WN stars while the broad-line stars, including 
IRS 13E2, could be H-free WN stars. This is consistent 
with He/H measurements ( [Najarro et al. 19'^ ) which cor¬ 
relate narrow-line stars with low He/H ( ^ 3) and bread¬ 
line stars with large He/H ( ^ 100). While LBVs generally 


stellar motion alone and thus would not appear as sym¬ 
metrical as, e.g., the Homunculus around rj Car. 

5.1. Lifetime of an LBV Nebula in the GC 

We can crudely estimate the “half-life” of an LBV shell 
(that is, how long before the two sides of the shell are 
differentially stretched by a factor of two if they are on 
separate circular Keplerian orbits) by 


t ~ 


2D 

Vi - Vo 


( 1 ) 


However, Paumard et al. (200l|) claim some or even all of 
these are actually not Hel stars and thus not likely WRs. 


where D is the diameter of the shell and Vi and Vo are the 
circular Keplerian velocity at the inner and outer edge of 
the shell, respectively. Assuming the object is near enough 
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to Sgr A* so that the gravitational potential of the super- 
massive black-hole is dominant ( ^ 3 pc; Genzel et al. 
1996|)7|the circular Keplerian velocity is: 


y^=100 km s-i , 


( 2 ) 


where R is in pc. Next, assuming D R, one finds that 
t is independent of D: 


t ~ 4 X 10"^yrs . 


(3) 


Thus, for IRS 13E2, any LBV shell is likely to exist for 
only ~ 2000 yrs after the star’s final molting. Since the 
circular orbital period of an object within a few parsecs of 
Sgr A* is 6 X 10"^yrs, the lifetime of any nebula will 
be less than an orbit, regardless of R. 

If the differential speed between the gas in the Bar and 
IRS 13E is a few tens of km s“^, a few thousand years is 
also the time-scale needed for IRS 13E to move from the 
ce ntre of the mini-cavity t o its present location (see Fig. 1 
in Coker fc Pittard 2000| ). Since the mm proper motion 
of the IRS 13E complex points back towards the mini¬ 
cavity ([Zhao fc Goss 1999 ) this raises the possibility of an 
association between the two: e.g., did IRS 13E carve out 
the mini-cavity as it traversed the Bar? 

Thus, there is strong evidence from the width of the 
He line at 2.058^m and the He/H ratio that IRS 13E2 has 


recentliy finished its LBV phase and is cuiienlly a biuad- 
line H-Iiee WN star. The pieseiice of LBV ejecta may also 
explain the high absorbing column and if true would place 
a strong constraint on the time elapsed since becoming a 
H-free WN star. 


6. An Extended Binary System 

In order for IRS 13E to be such a strong X-ray source 
while not being variable on time-scales of days or a year, 
it is likely to be a wide binary with a period significantly 
longer than a year. IRS 13E1, a possible massive dwarf O- 
star, is approximately 0.15" from 13E2 and is a possible 
companion since much of the mm emission lies between 


the separation between IRS 13E2 and IRS 13E1 corre¬ 
sponds to about 700 AU in projection, which, for a com¬ 
bined 70 Mq system, implies an orbital period of 2000 
years. At this separation the CWB shocks are likely to be 
largely adiabatic and would emit only weakly in the X-ray 
(Lx oc M'^v/d] [Stevens et al. 1992 ) unless the winds are 
particularly dense. 

However, on the plane of the sky, the IRS 13E complex 
is part of the Bar in the center of the mini-spiral. Thus, if 
the complex is located within or behind the dense gas of 
the Bar, it is possible that much of our estimated column 
is due to the Bar (which does not lie in front of Sgr A*) 
rather than the complex itself. In addition, the extinction 
towards the GC is known to be very patchy on small scales. 
This suggests that IRS 13E1 and IRS 13E2 may still be 
companions, but it would also imply that the LBV nebula 
has at least partially dispersed. 

The X-ray luminosities of CWB systems vary with 
time as well as depending on mass loss rates, wind ve¬ 
locities, and orbital separation. Models of WR 147, an 
extended CWB with a peri od of at least a few thousand 
years, predict Lx ^ O.ILq (Pittard et al. 2001), However, 
the mass-loss rate of IRS 13E2 is estimated to be more 
than ten times that of the WR 147 primary. Thus the 
X-ray luminosity of IRS 13E is consistent with a CWB. 

The proper motion of the mm sources in the IRS 13E 
complex have been measured to be ~ 250 km s“^ (Zhao & 
Goss 1998, 1999|). This is consistent with the interpreta- 


these two sources ( Zhao fc Goss 199g , 1999 ). IRS 13E3 is 
0.2" distant in projection and, although it has weak Pa- 
a/Hel emission so that it may be a WR-star in its own 
right, there is no mm emission directly between IRS 13E3 
and IRS 13E2. This implies that even if the three sources 
compose a triple system, IRS 13E3 is too far away to pro¬ 
duce significant emission due to colliding winds. 

A ver y crude es t imate of the intrin sic column in a 
CWB is ( Usov 1992 ; Stevens et al. 1992 ) 


Nr 


M 


~ 5 X 10 — cm 
dv 


(4) 


where the mass-loss rate M is in units of 10“^ Mq yr“^, 
the separation d is in units of 10^^ cm and the wind ve¬ 
locity V is in units of 1000 km s~^. For sensible values 
of M and v, if any significant fraction of our estimated 
column is intrinsic, d must be less than ^ 100 AU. But 


tion that they are high-velocity knots within surrounding 
stellar ejecta. Observations show that at least one of the 
mm sources has a spectral index suggestive of emission 
from ejected circumstellar material from a stellar enve¬ 
lope. 

In summary, the heavily absorbed X-ray spectrum of 
IRS 13E best matches a long period CWB system whose 
primary has recently exited its LBV phase. Long term 
monitoring at mm and K-band as well as long-integration 
Chandra observations will help determine precisely what 
type of binary is contained in the IRS 13E system. 
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